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ABSTRACT 

We present ID non-local thermodynamic equilibrium (non-LTE) time-dependent radiative- 
transfer simulations of a Chandrasekhar-mass delayed-detonation model which synthesizes 
0.51 Mq of ®®Ni, and confront our results to the Type la supernova (SN la) 2002bo over the 
first 100 days of its evolution. Assuming only homologous expansion, this same model repro¬ 
duces the bolometric and multi-band light curves, the secondary near-infrared (NIR) maxima, 
and the optical and NIR spectra. The chemical stratification of our model qualitatively agrees 
with previous inferences by Stehle et al., but reveals significant quantitative differences for 
both iron-group and intermediate-mass elements. We show that ±0.1 Mq (i.e., ±20 per cent) 
variations in ^®Ni mass have a modest impact on the bolometric and colour evolution of our 
model. One notable exception is the 17-band, where a larger abundance of iron-group elements 
results in less opaque ejecta through ionization effects, our model with more ^®Ni displaying 
a higher near-UV flux level. In the NIR range, such variations in ^®Ni mass affect the timing 
of the secondary maxima but not their magnitude, in agreement with observational results. 
Moreover, the variation in the I, J, and Ks magnitudes is less than 0.1 mag within ~10 
days from bolometric maximum, confirming the potential of NIR photometry of SNe la for 
cosmology. Overall, the delayed-detonation mechanism in single Chandrasekhar-mass white 
dwarf progenitors seems well suited for SN 2002bo and similar SNe la displaying a broad Si II 
6355 A line. Whatever multidimensional processes are at play during the explosion leading to 
these events, they must conspire to produce an ejecta comparable to our spherically-symmetric 
model. 
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1 INTRODUCTION 

The most widely-accepted model for Type la supe movae is the ther¬ 
monu clear disruption of a white dwarf (WD) star jHovle & Fowler! 
1 19601) in a binary system, although there is ongoing discussion 
about the combustion mode (pure deflagration or delayed detona¬ 
tion), the progenitor mass (Chandrasekhar mass or not), and the 
nature of the binary companion (another WD or a non-degenerate 
star). Ob servational evidence for diversity in the SN la popula¬ 
tion ( e.g., lvan Kerkwiik. Chang & Justhaml201(]l : lBadenes & Maoj 
[201^ seems to require multiple progenitor channels or explosion 
mechanisms. 

In a recent study, we found a promising agreement between a 
sequence of Chandrasekhar-mass delayed-de tonation models and 
observations of SNe la at maximum light iBlondin et alj [20l^ 
[hereafter B13]. The good match between our model DDC15 
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(~0.5 M© of ®®Ni) and the standard SN 2002bo motivates an in- 
depth study of this model at earlier and later times. _ 

Supernova 2002bo was first studied by iBenetti et alj )l2004l) . 
who favour a delayed-detonation explosion for this event based 
on the presence of intermediate-mass elements (IMEs) at high 
velocities and the lack of spectral signatures of CI/II associ¬ 
ated with unb umt carbon a t early times. Subsequent spectroscopic 
modeling by IStehle et'H] i2005h revealed a chemical stratifica- 
tion qualitatively similar to the fas t-deflagration W7 model of 
iNomoto. Thielemann & Yokoil l ll984l) . but with an offset to higher 
velocities and strong mixing. 

A striking feature of SN 2002bo is the large Doppler width of 
the Si II 6355 A line around m aximum light , which places this SN 
in the “broad-line” subclass of lSranch et alj l l2006l) n Such a broad 


^ Given the associated large blueshift at maximum absorption (> 
12000 kms“^ at maximum light), this SN is also part of the “high- 
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Si II line is systematically predicted in maximum-light spectra of 
the delayed-detonation models presented in B13. 

In the next section we hriefly present the numerical setup (hy¬ 
drodynamics and radiative transfer). The observational data and 
our method of analysis are presented in Section [3 We then con¬ 
front the predicted bolometric (SectionlUl, colour (Section]^, and 
spectroscopic (Section|^ evolution of our model to observations of 
SN 2002bo, both in the optical and NIR. We illustrate the sensitiv¬ 
ity of our model results to ±0.1 Mq (i.e., ±20 per cent) variations 
in ®®Ni mass in Section|7] We compare the abundance distributions 
of various e lements in o u r inpu t hydrodynamical model with those 
inferred by IStehle et'H] ( l2005l) in Section]^ A discussion on the 
nature of broad-lined SNe la and conclusions follow in Section|^ 


2 NUMERICAL SETUP 

The hydrodynamics of the explosion and the rad iative-transfer 
treatme nt is analogous to that presented in B13 and lPessart et alj 
( l2014ch [hereafter D 14c]. We refer the reader to those papers for an 
in-depth discussion of our numerical setup. As with all our previous 
SN studies, our models are iterated until convergence is achieved at 
all depths in the ejecta (see Appendix I aIi. 

We use the Chandrasekhar-mass delayed-detonation model 
DDC15 presented in B13, and summarize its basic properties 
in Table [T] Unlike B13, we apply a small radial mixing to 
the hydrodynamical input with a characteristic velocity width 
Aumix = 400 kms~^ to s mooth sharp variations in composi¬ 
tion (see lPessart et alj|2014al : hereafter P14a). Moreover, we con¬ 
sider additional radioactive decay chains (only ®®Ni decay was 
considered in B13), resulting in non-negligible energy input be¬ 
yond 'xjISOOO kms“^ where the ejecta is deficient in ®®Ni (see 
model PPC10_A4P1 in P14a). More importantly, including the 
'*®Cr —>■ —>• "‘^Ti decay chain leads to a ~ 100-fold increase 

in the Ti mass fraction in the inner ejecta {v < 13000 kms“^; 
see Fig.O, which influences sizeably the opacity in the B-band at 
maximum light and beyond. 

The chemical stratification is typical of standard delayed- 
detonation models. Iron-group elements (IGEs; Sc to Ni, Mtot ~ 
0.7 M 0 ) dominate below ~10000 kms“^. Pespite the presence 
of ®®Ni at <2000 kms~^ due to our imposed radial mixing, the 
innermost ejecta are dominated by stable IGEs. Intermediate-mass 
elements (Na to Ca, Mtot ~ 0.6 Mq) are most abundant between 
^10000 kms“^ and ~30000 kms“^. The outermost ~ 0.01 Mq 
of the ejecta is dominated by the composition of the progenitor WP 
star (C/0, with traces of ^^Ne, and solar abundances for all other 
elements). 

The ejecta density profile at the start of our radiative-transfer 
calculations (0.5 days past explosion) is shown in Fig.[2 At these 
times, UV photons still interact with the ejecta material in shells 
moving at a velocity > 40000 kms“^, so we linearly extrap¬ 
olate the hydrodynamical input beyond ~45000 kms“^ up to 
~70000 kms“^, as done in P14a (dotted line in Fig. O. The 
density profile is reminiscent of standard delayed- detonation mod¬ 
els (see, e.g.. iKhokhlov, Muller & Hoflichl Il993h . with a quasi¬ 
exponential decrease with velocity. Also shown is the density pro¬ 
file of the W7 model used as a reference in many SN la studies. It 
is nearly indistinguishable from our model below ~ 10000 km 

velocity” s ubclass of|Wang_gl_alJ l2009h and the “high-velocity-gradient” 
subclass of iBenetti et alj ]2005h . 


Lagrangian Mass [Mq] 
0.1 0.5 1.0 1.3 1.4 





Velocity [10^ km s 

Figure 1. Abundance profiles in velocity space for C/O/Ne (top), and se¬ 
lected IMEs (middle), and IGEs (bottom; the profiles for ®®Ni and '^®Cr 
at this time are shown as a dashed and dotted line, respectively) in model 
DDC15 at the start of our radiative-transfer simulations (0.5 days past ex¬ 
plosion). The upper abscissa shows the Lagrangian mass coordinate. 


(or within a Lagrangian mass coordinate of ~1 M©). The key dif¬ 
ference lies in the amount of mass contained at high velocities: It 
is a factor of ~2 larger beyond -^15000 kms“^ in PPC15, re¬ 
sulting in a ~20 per cent larger kinetic energy (-^ 1.5x10®^ erg 
cf. ~1.2x 10®^ erg for W7; see Table Q}. The composition is also 
drastically different, with IMEs confined to < 15000 km s“ ^, while 
they are abundant beyond 20000 kms“^ in PPG 15. 

The radial mixing applied to the hydrodynamical input and 
the extended radial grid result in small differences in abundance 
between the PPG 15 model presented in B13 and the version used 
in this paper. In particular, the initial ®®Ni mass for this model is 
now ~ 0.51 Mq, cf. ~ 0.56 Mq in B13 (see their Table 1). These 
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Lagrangian Mass [Mq] 

0.10.5 1.0 1.3 1.38 W7 


0.1 0.5 1.0 1.3 1.38 1.4 DDC15 



Figure 2. Density profile at the start of our radiative transfer simulations 
(0 .5 days past explosion; black line), com pared to the W7 model (grey line) 
of iNomoto. Thielemann & Yokoil jl984) . The solid dots conespond to the 
initial radial grid used in CMFGEN. The dotted curve shows the extrapolated 
density profile beyond ~45000 km (see text). The upper axis gives the 
Lagrangian mass coordinate for both DDC15 (black) and W7 (grey). 


differences only have a modest impact at bolometric maximum: A 
bolometric rise time of 17.57 d (cf. 17.76 d in B13), a peak bolo¬ 
metric luminosity of 1.14x10'^® ergs“^ (cf. 1.22x10'^® ergs“^), 
and a.B — R colour of 0.04 mag (cf. —0.02 mag). 

In what follows, we use model DDC15 to discuss SN 2002bo 
and broad-lined SNe la in general. However, we also study the sen¬ 
sitivity of our results by comparing the predicted light curves and 
spectra of delayed-detonation models endowed initially with ±20 
per cent of mass of ®®Ni (models DDCIO and DDC17 in Table[T]l. 
This comparison is particularly instructive to understand the degen¬ 
eracy (or lack thereof) of SN la properties. 


3 OBSERVATIONAL DATA AND METHODOLOGY 


Our d ata on SN 2002bo are mostly take n from Bgnettj_et_alJ 
(l2004l). with additional opti cal spectra from iBlondin et alj 1 2012h 
and Isilverman et ^ ( l2012h. We use the near-i nfrared (JHKs) 
photometry published by Krisciunas et alj ( |2004|) . We select spec¬ 
tra whose relative flux calibration agrees with photometry-based 
colours to within 0.1 mag. As in B13, we adopt a recession ve¬ 
locity cz = 1293 kms“^ for the host galaxy, a distance modu¬ 
lus fi = 31.90 ± 0.20 mag and a total reddening E{B — V) = 
0.41 ± 0.07 mag. 

The method for generating pseudo-bolometric light curves 
for both models and data is analogous to that presented in B13. 
We specify each time the bluest and the reddest band used in the 
integration, such that Lu^k, corresponds to the luminosity (in 
erg s“^) obtained by integrating UBVRIJHKs magnitudes. 

Since there is no good definition of a photosphere through 
opacity means in SNe la, we characterize the spectrum-formation 
region using the spatial distribution of the optical (3000-10000 A) 
flux instead. The fractional contribution to the total flux at wave¬ 
length A from a given ejecta shell at radius r and thickness dr at 
time t since explosion is: 



Figure 3. Illustration of the wavelength dependence of the spectrum- 
formation region, represented here using the velocity above which 50% of 
the flux at a given wavelength is emitted. This velocity is computed based 
on equation but replacing the wavelength-averaged flux contribution 
SF{v) of equation ^2) with the non-averaged SF{r, A, t) of equation (T). 


SF{r,X,t) = ^ j Azr]{p,z,X,t) e p dp, (1) 

where D is the physical distance, Az is the projected shell thick¬ 
ness for a ray with impact parameter p, p is the emissivity along the 
ray at p and z, and r is the ray optical depth (at A) at the ejecta loca¬ 
tion (p, z). Because of homologous expansion (v = rjt), surfaces 
of constant projected velocity are planes perpendicular to the 2 : di¬ 
rection. Hence, the quantity SF{r, X, t) mixes contributions from 
regions with the same radial velocity but with very different pro¬ 
jected velocities, hence Doppler shift. The averaged contribution 
of a shell with velocity v (corresponding to r ± dr 12) over the 
wavelength interval [Ai, A 2 ] is given by: 


5F{v) = 


J^^SF{r,X,t)dX 

|Ai — A2I 


( 2 ) 


A representative velocity for the spectrum-formation region in the 
wavelength range Ai to A 2 may be the velocity above which 50% 
of the flux in that range is emitted, noted V 1 / 2 , such that: 


E SF{v) 

v^Vl/2 _ 1 

_ _ r\ 5 

E SF{v) 2 

all V 

where the sum in the denominator runs over the entire velocity grid 
spanned by our simulation. 

In what follows this velocity is only used for illustration pur¬ 
poses, as the averaging in equation id occults the strong variation 
of v-i /2 with wavelength (Fig.d- This variation reflects the dom¬ 
inanc e of line opacity in SN la ejecta (see, e.g.. |Pinto & EastmanI 

l200flh . 


4 BOLOMETRIC EVOLUTION 

The pseudo-bolometric light curves for both our DDC15 model 
and SN 2002bo are shown in Fig. |4] Our model predicts a rise to 
bolometric maximum of ~17.6 d, while the first observations of 
SN 2002bo are ~14 days before peak. At these early times, the 
bolometric evolution is well matched by our model, albeit with 
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Table 1. Model parameters and nucleosynthetic yields for selected species at the stait of our radiative-transfer calculations (0.5 d past explosion). The 
mass is given at texp ~ 0. Models DDCIO and DDC17 will be discussed in Section]?] 


Model 

ptr 

£)kin 

t)(56Ni) 

®®Nit=o 

Ni 

Co 

Fe 

Ti 

Sc 

Ca 

S 

Si 

Mg 

Na 

0 

C 


[g cm-3] 

[ergs-1] 

[kms-l] 

[Mq] 

[Mq] 

[Mq] 

[Mq] 

[Mq] 

[Mq] 

[Mq] 

[Mq] [Mq] 

[Mq] 

[Mq] 

[Mq] 

[Mq] 

DDC15 

1.8(7) 

1.465(51) 

1.12(4) 

0.511 

0.516 

3.44(-2) 

0.114 

O 

00 

4.53(-2) 

0.197 

0.306 

1.14(-2) 

1.68(-5) 

0.105 

2.73(-3) 

DDCIO 

2.3(7) 

1.520(51) 

1.16(4) 

0.623 

0.622 

4.11(-2) 

0.115 

1.10(-4) 

2.90(-8) 

4.10(-2) 

0.166 

0.257 

9.95(-3) 

1.25(-5) 

0.101 

2.16(-3) 

DDC17 

1.6(7) 

1.459(51) 

1.08(4) 

0.412 

0.421 

2.84(-2) 

0.112 

1.14(-4) 

6.16(-8) 

4.73(-2) 

0.222 

0.353 

1.79(-2) 

2.47(-5) 

0.152 

3.80(-3) 


Notes: Numbers in pai'enthesis con'espond to powers of ten. The deflagration velocity is set to 3% of the local sound speed ahead of the flame for all models; 
ptr is the transition density at which the deflagration is artificially turned into a detonation; £’kin is the asymptotic kinetic energy; i;(^®Ni) is the velocity of 
the ejecta shell that bounds 99% of the total ^®Ni mass. 


Days since Explosion (DDC15) 



Days from Bolometric {U Kg) Maximum 


Figure 4. Top: Pseudo-bolometric {U —>■ Ks) light cuiwe for model DDC15 (black) and SN 2002bo (open circles). The error bars take into account measure¬ 
ment, extinction, and distance errors. The grey dashed line shows the instantaneous rate of decay energy. Coloured lines correspond to models with ±0.1 Mq 
of ^®Ni compared to DDC15. We show close-up views of the early-time B ^ I (left inset) and late-time U ^ I (right inset) pseudo-bolometric light curves, 
due to the lack of [/-band or NIR data at these times. The upper abscissa gives the time since explosion for model DDC15. Bottom: Contribution of the 
U Ks (dash-dotted line), U —y I (dashed line) and B I (dotted line) pseudo-bolometric luminosities to the total bolometric luminosity. Also shown 
are the contributions of the UV (A < 3000 A; grey solid line) and IR (A > 1 p-m; grey dashed line) ranges. 
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Table 2. Bolometric properties. Models DDCIO and DDC17 will be dis¬ 
cussed in Sectionl?] 


Model 

56Ni 

imq] 

^rise 

[d] 

OL 

-^peak 

[ergs“l] 

AMi5 

[mag] 

Q 7 

•^dep 

DDC15 

0.511 

17.57 

3.4 

1.14(43) 

0.72 

1.04 

0.95 

DDCIO 

0.623 

17.13 

3.3 

1.38(43) 

0.74 

1.00 

0.95 

DDC17 

0.412 

18.60 

3.3 

9.10(42) 

0.67 

0.99 

0.97 


Notes: Numbers in parenthesis correspond to powers of ten; a is the 
exponent of the power-law fit to the early-time (texp 3 d) luminosity, 
L{t) oc AMi 5 is the bolometric magnitude decline between 
bolometric maximum and 15 d after; Q-, is the ratio of the peak bolometric 
luminosity (ipeak) to the instantaneous decay luminosity; Tyjep the 
fraction of decay energy at peak actually deposited in the ejecta. 


a slightly larger luminosity (but still within the errors). This sug¬ 
gests the ®®Ni distribution is adequate, in particular its presence 
at >15000 kms“^. In W7, the ®®Ni mass fraction drops to 10“® 
at ~ 12000 kms“^, and the early-time lum inosity is significa ntly 
lower than observed in SN 2002bo (see also lStehle et alfcoOSh . 

A power-law fit L{t) ex to the early (fexp ^ 3 d) bolo¬ 
metric light curve of model DDC15 yields a = 3.4 ( Table [Hi. De¬ 
viatio ns from a “fireball” model (for which a = 2; [Nugent et alj 
[201 ih are expected theoreti cally and confirme d observationally, 
as in the recent SN 2013dy jZheng et alj[2013h . Variations in the 
early-time luminosi ty can result from outward mixing of ®®Ni (e.g., 
[Piro & Naka3[20l4 . or from structural changes in the outer ejecta, 
as in the pulsational-delayed-detonation models of D 14a. 

The lower panel of Fig. [4| gives the fractional contribution to 
the total bolometric luminosity of various filter sets, as well as in 
the UV (< 3000 A) and IR (> 1 /rm) ranges. At maximum light, 
we recover 94 per cent of the true bolometric luminosity with our 
U —>■ Ka pseudo-bolometric light curve (see also Table [IB. The 
remaining 6 per cent is mostly emitted blueward of the (7-band. 

From bolometric maximum until the onset of the radioactive 
tail -^35 days later the 7 -ray escape fraction increases from to 
-^50 per cent. The decline in bolometric magnitude between maxi¬ 
mum light and 15 days later (noted AM 15 in Table[B is 0.72 mag, 
comparable to that inferred for SN 2002bo based on its U —>■ Ks 
pseudo-bolometric light curve (^0.85 mag). 

Past maximum, the transition to an optically-thin ejecta is con¬ 
comitant with a shift in ionization and a greater release of radi¬ 
ation energy, causing a small bump in the bolometric light curve 
around ~20 days past maximum light (~30 days past maximum in 
SN 2002bo; see Fig.0. 

Past -^30 days from maximum light (~50 days since explo¬ 
sion), the bolometric luminosity follows the instantaneous rate of 
energy deposition. The 7 -ray escape fraction continues to increase, 
reaching ^80 per cent at ~100 days past explosion. The late¬ 
time (> 50 days past maximum) pseudo-bolometric (U —>■ 7) 
decline is compatible with SN 2002bo (-^0.024 mag day cf. 
~0.022±0.007 mag day-1). 

The kinetic energy, ®®Ni mass, and ejecta mass of our model 
yield a bolometric light curve consistent with observations of 
SN 2002bo at all times. This represents a compelling check for 
the standard Chandrasekhar-mass delayed-detonation model for 
this class of SN la light-curve morphology. Moreover, the good 
match with a ID model suggests that asymmetries, if present, may 
only have a minor role in shaping the bolometric light curve of 
SN 2002bo. 
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Figure 5. Evolution of the gas temperature profile (top) and of the cobalt 
ionization fraction (bottom) in model DDC15, between 1 and 100 days 
since explosion (note the logarithmic time steps). Solid dots mark the ve¬ 
locity location above which 50% of the optical flux is emitted, which we 
use to characterize the spectrum-formation region (see Section[B. 


5 PHOTOMETRIC EVOLUTION 

The colour evolution of SNe la reflects both the global evolu¬ 
tion of ejecta properties (temperature, mean ionization) and the 
wavelength-dependent opacity, and is thus very difficult to accu¬ 
rately model. 

Initially, all bands experience a steep brightening due to the 
increase in radius. The brightening is steeper for the bluer bands 
(see the a exponent in Table[B due to the temperature increase in 
the spectrum-formation layers (solid dots in Fig. [B, which shifts 
the optical and NIR colours shown in Fig.[Bto the blue. Because of 
the lower opacity/emissivity compared to the optical, the NIR flux 
starts to decrease as the spectral energy distribution (SED) hardens, 
as early as 7 days prior to bolometric maximum for the H and Ka 
bands in model DDC15 (Fig.0. In SN 2002bo, the initial peak in 
these bands occurs ~3 days later than in our models, similar to the 
evolution in the J band (which is well matched by our model). 

In the optical, the shorter rise times in the U and B bands with 
respect to V and R reflects the temperature decline in the spectrum- 
formation region after about one week past explosion. By maxi¬ 
mum light, most of the flux originates from layers <15000 kms”^ 
where the temperature is decreasing (Fig.B- 

The U — V and B — V colours shift to the red at ~ 10 days past 
explosion (one week before B-band maximum), i.e. 10-20 days 
before the other colours shown in Fig.[B This earlier shift in colour 
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Days Since Explosion (DDC15) 

0 20 40 60 80 100 20 40 60 80 100 20 40 60 80 100 20 40 60 80 100 

-20 

-19 

-18 

-17 

-16 

-15 

-20 

-19 

-18 

-17 

-16 

-15 

0 20 40 60 80 0 20 40 60 80 0 20 40 60 80 0 20 40 60 80 

Days from Bolometric (!/—>■ Ks) Maximum 

Figure 7. Optical (UBVRl) and NIR (JHKs) light curves for model DDC15 (solid lines) compai'ed to SN 2002bo (dereddened; open circles). The eiTor 
bars take into account measurement, extinction, and distance eri'ors. Coloured lines correspond to models with diO.l Mq of ^®Ni compared to DDC15. The 
upper abscissa gives the time since explosion for model DDC15. 



Table 3. Light-curve properties. For each parameter, the upper line corresponds to model DDC15 (the upper and lower values, when added to the DDC15 
value, correspond to models DDCIO and DDC17, respectively), and the lower line to SN 2002bo (Icr en'ors are given in between parentheses). 


Parameter 

U 

B 

V 

R 

I 

J 

H 

Ks 

o in L(t) cc 
(t^exp < 3 d) 

A n+ 1-2 

^■^+ 1.4 

5 Q“0.1 

^•^- 0.3 

trr Q —0.6 
O.O_ 0.4 

4 1 - 0 -® 
^■^- 0.3 

P 7-0.5 
'^•'-0.2 

2-6+1 

2 

^•'^-0.1 

2 5—00 
'^•^-0.0 

^^peak 

[day] 

-1 in-o -65 
- 0.68 ( 0 . 70 ) 

U.oy_ 0.65 
+ 0.61 ( 0 . 60 ) 

_ 1_3 

-|-o.za_ 0_07 

+ 2.44 ( 0 . 67 ) 

_ 1_2 1 q+ 0.39 
-r^-J-y_o .99 

+ 2.15 ( 0 . 54 ) 

,1 rri-0.53 
-|-1.00_Q 29 

- 1.47 ( 0 . 34 ) 

- 4 . 05 +: 6 l 
-1.39 ( 0 . 14 ) 

0-00^1 30 

- 1 . 10 ( 0 . 14 ) 

—6 62 “*^'®® 
D.oz^_o 68 

- 0.84 ( 0 . 28 ) 

max 

[day] 





+ 20 . 4211 " 
+ 24.45 ( 1 . 52 ) 

+ 23 . 421+3 
+ 28.36 ( 0 . 39 ) 

+2O.66I3 g3 
+ 24.93 ( 0 . 42 ) 

+ 20 . 84 l 3 ;|l 
+ 24.81 ( 0 . 98 ) 

-^peak 

[mag] 

iy. 4 U^ 0.59 
- 19.38 ( 0 . 22 ) 

— 10 1 
iy.j.o,|_0_22 

- 19.20 (O.ll) 

- 19.20 ( 0 . 10 ) 

iy.zu,^ 0 .i 3 

- 19 . 10 ( 0 . 09 ) 

- 18 . 83 ;°:“^ 
- 18.75 ( 0 . 10 ) 

18 c;q+0-02 

±o.oo_Q Q3 

- 18.43 ( 0 . 14 ) 

1 Q 1 q+0.02 

10 . 13 _o.o 7 
- 18.18 ( 0 . 14 ) 

- 18 . 071 +° 

- 18 . 13 ( 0 . 14 ) 

'^2naniax 

[mag] 





-18-861^0 
- 18.33 ( 0 . 13 ) 

_17 03+0.14 
± 1 •»a_o Q2 

- 17.44 ( 0 . 13 ) 

1 ft Qp+ 0.06 

10 . 03 ^_o.oi 
- 18.21 ( 0 . 14 ) 

1 ft 7Q + O.O6 

10 ./ 3 ^ 0.01 
- 18.05 ( 0 . 15 ) 

AMi5 

[mag] 

1.53 ( 0 . 13 ) 

-L-DO^O .24 

1 . 18 ( 0 . 07 ) 

n 71 — 0.05 
'^•'-‘-+ 0.17 

0.70 ( 0 . 03 ) 

n 54+0 08 

U. 04 _o 03 

0.71 ( 0 . 04 ) 

0 . 08 l°;?^ 

0.64 ( 0 . 03 ) 

1 Ql —012 
-‘-•'^-‘--0.11 
1.65 ( 0 . 06 ) 

o.oein^ 

0.37 ( 0 . 03 ) 

tJ.U 4 _ 0.06 

0.42 ( 0 . 07 ) 

AM/At 5 o+ 
[mag day ”1] 

U. 02 ^ 4 ,^ 0.000 
0 . 019 ( 0 . 013 ) 

n ni 

U.UID^O.OOO 
0.011 ( 0 . 005 ) 

U.UZD_o 001 
0.023 ( 0 . 003 ) 

n QQn+O OOl 

U.UCiU_0.001 
0.033 ( 0 . 002 ) 

0 nQ4+0-002 

U.U04_q QQ2 
0.041 ( 0 . 001 ) 

U.U 0 D_Q QQ2 

n 041 

U.U 41 ^ 0.000 

0 0*37—0-002 

u.uo/^_o.oo2 


Notes: the /-band light curves for models DDCIO and DDC15 have several glitches around maximum light, hence the parameter values ai'e less reliable; 
Afpeak is the time difference between maximum light in a given band (first maximum for IJHKs) and pseudo-bolometric (U Ks) maximum; 
Af 2 ndmax i^ ^i^^ difference between the secondary maximum (IJHKs only) and pseudo-bolometric (U Ks) maximum; AM/Af 5 o+ is the 
average magnitude decline rate past -1-50 d from maximum light. 
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Days since Explosion (DDC15) 

0 20 40 60 80 100 



Days from Bolometric {U Kg) Maximum 

Figure 6. Comparison of colour curves for model DDC15 (black) and 
SN 2002bo (dereddened; open circles). The en'or bars take into account 
measurement and extinction errors. Coloured lines correspond to models 
with ±0.1 Mq of ®®Ni compared to DDC15. Note that the redder colour 
curve corresponds to the model with less ®®Ni (DDC17). The upper ab¬ 
scissa gives the time since explosion for model DDC15. 


results from the increase in opacity from once-ionized IGEs that 
deplete the flux in U and B relative to redder bands. 

The Ami 5 {B) parameter for our model is 1.03 mag, con¬ 
sistent with that inferred for SN 2002bo (1.18 ± 0.07 mag; see 
Table O. In the redder optical bands and in the NIR range, the 
post-maximum decline rate is sensitive to the location of line emis¬ 
sion. In the IHKg bands, line emission is sufficient to maintain 
a near-constant flux level despite the overall decline in bolomet¬ 
ric luminosity HKs) « 0 mag), while the weaker line 


emission within the J band causes a more pronounced decline 
(AMi 5 (J) « 1.3 mag). 

At late times, the colours shown in Fig. [6] become progres¬ 
sively bluer despite the near-constant temperature (7000-8000 K 
in the ejecta layers between 5000 kms“^ and 15000 kms“^; see 
Fig. [5): The low ejecta densities favour the formation of forbid¬ 
den lines, whose wavelength location impacts the colour evolution 
independent of temperature. 

Our model reproduces the complex evolution of the optical 
and NIR colours of SN 2002bo out to ~80 days past maximum 
light. This gives us confidence that the apparent red colours of 
SN 2002bo are indeed due to extinction in the host galaxy. More¬ 
over, the ob served colour curves are corrected for reddening using 
the classical fCardelli, Clayton & MathisI l ll989l) extinction law with 
a standard value of the total-to-selective extinction ratio Rv = 3.1. 

However, the model is very red within the first ~5 days 
since explosion, where we lack SN 2002bo d ata. Earlier observa- 
tions of sev eral other SNe la ( e.g., SN 201 Ife. lNugent et aT]|201ll ; 
SN 2013dv. lzheng et al.l2013l) . all of which display a nan'ower Si II 
6355 A line than SN 2002bo around maximum, reveal significantly 
bluer colours than predicted by our model. Such blue colours at 
early times are difficult to reproduce with standard ID delayed- 
detonation models, even with strong outward mixing of ®®Ni, and 
instead appear to require some form of hydrodynamical interac¬ 
tion to heat up the outer ejecta layers, as in the pulsating delayed- 
detonation models of D14a0 Earlier observations of a broad-lined 
SN la such as SN 2002bo would enable to check the prediction of 
red early-time optical colours in our model. 


6 SPECTROSCOPIC EVOLUTION 

We describe below the evolution of the optical (Fig. (8)1 and NIR 
(Fig.|9) spectra of our model and compare them to SN 2002bo ob¬ 
servations. Tables with detailed line identifications, and a collection 
of plots illustrating the influence of individual ions on the synthetic 
spectra, are available in AppendiceslOandlDl respectively. 

6.1 Early evolution 

During the first -^10 days past explosion, the majority of the 
optical flux emerges from the IME-rich regions of the ejecta at 
>15000 kms“^. The optical spectrum is thus dominated by broad 
lines of Mg II, Sill, SlI, and Call. We also note the presence 
of 0 17773 A, which causes a deep absorption blueward of the 
Ca II 8500 A triplet. Iron-group elements, most notably Sc II, Ti II, 
Cr II, Fe II, Co II, and Ni II, also contribute to the flux shortward of 
-5000 A. 

We obtain high-velocity absorption (> 25000 kms“^ 

blueshifts) in both the Call H & K and 8500 A triplet features. 
These remain optically thick well beyond 30000 km where the 
Ca abundance is at its primordial solar value (X « 6 x 10“®), 
but their blue wings are affected by absorption from other lines and 
thus cannot easily serv e as probes of the progenitor metallicity (see 
also lUentz et al.|[2000l in the context of the W7 model). The Sill 
3858 A line has negligible impact on our synthetic spectra, and 
the predicted double-absorption feature at —3500 A seen, e.g., in 
the — 7.7 d spectrum results from Call alone. Multiple absorption 

^ Although such an interaction serves to decelerate these outer ejecta lay¬ 
ers, inhibiting the formation of broad spectral lines as in SN 2002bo. 
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Figure 8 . Optical spectroscopic evolution of model DDC15 (red) compared to SN 2002bo (black), between —12.9 d and +83.6 d from pseudo-bolometric 
{U —>■ Ks) maximum. The tickmai'ks on the ordinate give the zero-flux level. The observed spectra have been de-redshifted, de-reddened, and scaled to match 
the absolute V-band magnitude inferred from the corresponding photometry. An additional scaling (F’scT green label) has been applied to the synthetic spectra 
to reproduce the mean observed flux in the range 5000-6500 A. 
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components in both Ca II features are routine ly observed in early- 
time SN la spectra fe.g. JChildress et a02014h . 

Unburnt material is located in the outermost 
(>20000 kms“^) low-density layers of the ejecta. This con¬ 
figuration makes the associated spectral signatures of carbon (e.g., 
ClI 6578 A) too weak to be visible. The Ol 7773 A line can in 
principle serve to constrain the amount of unburnt oxygen, but it 
forms below ^25000 kms“^ in our model and hence results in 
part from C-buming ash. 

In the NIR range, our model reproduces the overall shape of 
the earliest spectrum of SN 2002bo at —8.5 d (~9 days past explo¬ 
sion based on model DDC15), but predicts a higher flux level than 
observed (Fsd = 0.64). Three ions are responsible for the spec¬ 
tral signatures seen: O I, Mg II, and Si II (see Fig. C2, upper-right 
panel). The two most notable features at these early times are a P- 
Cygni profile due to the Mg II 10927 A triplet and a broad emission 
feature at ~ 1.7 fim due to the Si II 16931 A doublet. 

Lines from IGEs are completely absent from the NIR range 
at these early times, although lines of Sc II, Ti II, Cr II, and Fe II- 
III (with weaker absorptions from Coll and Nill in the (7-band) 
contribute to the optical flux. The majority of the NIR flux emerges 
from ejecta layers <20000 kms“^ where most IGEs are at least 
twice-ionized, and contribute opacity at much shorter wavelengths. 

The earliest spectrum at ~13 days before bolometric maxi¬ 
mum, however, shows that we predict too little flux around 4000 A, 
due to excess absorption from Sc II 4247 A and several Ti II lines, 
and that the blue extent of certain features is not well matched, in 
particular the Ee II absorption complex around 4800 A, and the two 
Si II 5972 A and 6355 A features. The 01 7773 A absorption is also 
too strong in our model. These discrepancies are related to the low 
ejecta ionization fraction at 20000 km (Eig.[5j, and suggest a 
stronger outward mixing of ®®Ni than adopted here. Another pos¬ 
sibility is that our extrapolated density profile at large velocity is 
steeper in reality (see Fig.|^. 

6.2 Maximum-light phase 

Around maximum light the optical spectra are still dominated by 
strong lines of singly-ionized IMEs, with IGEs absorbing the flux 
at <5000 A. The drop in the Mg and O mass fractions below 
~ 15000 kms“^ results in a gradual weakening of the Mg II and O I 
7773 A lines and reveals the transition from the C-buming to the 
0-burning regime. We note the disappearance of the high-velocity 
absorption associated with Ca II H&K and the 8500 A triplet, as 
these transitions turn thin beyond ^^20000 kms“^. 

In the near-infrared, we note the appearance of emission fea¬ 
tures associated with IGEs, most notably permitted Co II and for¬ 
bidden Com transitions. There are no published NIR spectra of 
SN 2002bo around maximum light to compare against our model. 

6.3 Post-maximum and transition to nebular phase 

From ~10 until <40 days past maximum light the majority of 
the flux emerges from the outer IGE-rich regions of the ejecta 
(<10000 kms“^). Some strong lines of IMEs survive in the optical 
range (Si II, S II, Ca II), although they become weaker and narrower 
with time. Most of the spectrum is now shaped by lines of IGEs 
such as Ti II, Cr II, Fe II-III, and Co II-III. 

The SED becomes redder due to the gradual drop in tempera¬ 
ture in the inner ejecta layers. Radiation cooling dominates over ex¬ 
pansion cooling, with a growing importance played by forbidden- 
line transitions. A notable example is the 5900 A feature, well 



Wavelength [A] 

Figure 9. Near-infrared spectroscopic evolution of model DDC15 (red) 
compai'ed to SN 2002bo (black), between —8.5 d and +84.6 d from pseudo- 
bolometric (U —>■ Ks) maximum. The flux has been scaled by for bet¬ 
ter visibility, and the tickmarks on the ordinate give the zero-flux level. The 
observed spectra have been de-redshifted, de-reddened, and scaled to match 
the absolute J-band magnitude inferred from the corresponding photome¬ 
try, when available. An additional scaling (-Fgcli gi'een label) has been ap¬ 
plied to the synthetic spectra to reproduce the mean observed flux in the 
range 10000-11500 A. 

matched by our model, and associated with [Co III] 5888 A. This 
feature is often mistaken for Na I D (see discussion in lPessart et all 
l2014hh . 

This cooling further enhances absorption by Cr Il/Ti II transi¬ 
tions around 4000-4300 A. Several Cm lines are responsible for 
the rapid change in relative flux around ~5300 A between - 1 - 16.2 d 
and 5-28.2 d (Fig. [U, which reflects the rapid evolution of the 
Cr II/III ratio. 

In the NIR, IIWII recombination of IGEs causes a gradual 
weakening of Co III lines (primarily forbidden transitions) and a 
strengthening of Coll lines (primarily permitted). In some bands 
(H and Ka), emission by Coll is sufficient to maintain a near¬ 
constant flux level and causes an earlier shift to the red of the 
V — Ka colour curve with respect toV — J (Fig. [Hi. 

The strengthening of the Co II emission is partly responsible 
for the mismatch in the Ca II 8500 A triplet emission profile be¬ 
tween our model and SN 2002bo. Co II lines blueward of ~8500 A 
emit photons in the inner ejecta (<10000 kms“^) that are red- 
shifted into resonance with the Call 8500 A triplet transitions in 
higher-velocity shells (see also D14c). The resulting boost in emis- 
sivity explains the slower decline in the /-band in our model com¬ 
pared to SN 2002bo (cf. Fig.|7]and Table|3j. 
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6.4 Late-time evolution out to ~100 days past explosion 

From ~40 days past maximum onwards the spectrum-formation 
region is almost exclusively located in the IGE-rich layers of the 
ejecta helow 10000 km s“^. Apart from the two strong Call H&K 
and 8500 A triplet features and traces of Si II 6355 A, which retain 
their P-Cygni profile morphology, the only optical lines from IMEs 
are optically-thin transitions of [Sll], [Sill], and [Arm]. Beyond 
l/im, the only emission from IMEs comes from [S II] and Call 
transitions, although the impact on the observed flux is negligible 
given the overlap with stronger emission from Ee II and Co II. 

Non-thermal ionization prevents a complete III—>II recombi¬ 
nation of IGEs in the inner ejecta, resulting in increasing emis¬ 
sion from lines of [Eelll] and [Com], now favoured due to the 
low ejecta density. Some of these lines represent a fair fraction of 
the total flux in a given band: At -^100 days past explosion, the 
[Felll] 4658 A and [Co III] 5888 A lines alone contribute ~15 per 
cent of the total flux in the B and R bands, respectively. The ratio 
of the total flux emitted in both lines follows the increase in the 
Ee/Co abundance ra tio resulting from ®®Co —>■ ®®Ee decay (see also 
iKuchner et aljl99^ . The late-time colour evolution is thus largely 
determined by the wavelength location of individual strong transi¬ 
tions and no longer reflects the evolution of the ejecta temperature, 
which remains almost constant in the spectrum-formation region at 
these times (Eig.|5}. 

In the NIR, we predict the gradual emergence of [Nill] 
1.94 ^m, due to stable Ni in the innermost ejecta (< 2000 kms“^; 
see Fig. [T]l, resulting from burning at high densities characteristic 
of a C handrasekhar-mass WD progenitor (see, e.g.. iHoflich et alj 
1 20021) . However, most of the emission just blueward of 2 /rm 
is in fact due to lines of [Com] in our model, not [Ni II] 
as recently suggested by iFriesen eT iD i l2014h for SNe 2001fe, 
2 003du, and 2014J based on the ID del ayed-detonation models 
of lDomlnguez, Hoflich & Stranierd ( l200lh . Unfortunately, the NIR 
spectrum of SN 2002bo at -1-56.3 d does not cover this feature. 

The match between the spectra of model DDC15 and 
SN 2002bo over the first -^100 days of its evolution suggests the 
composition of the ejecta in our model is adequate at all depths. It 
also shows our time-dependent radiative-transfer treatment consis¬ 
tently predicts the evolution of the gas properties, in particular its 
ionization state. 


7 SENSITIVITY TO Ni MASS 

Our model DDC15 seems well suited to match the observed proper¬ 
ties of SN 2002bo. In this section, we investigate whether delayed- 
detonation models with ± 0.1 M© (i.e., ±20 per cent) variations in 
®®Ni mass also provide a satisfactory match to this event (models 
DDCIO and DDC17 of B13; see TableO. This is of particular rel¬ 
evance to assess the actual accuracy of abundance determinations 
in SNe la and to understand the degeneracy (or lack thereof) of 
their properties. We discuss in turn the impact on the early bright¬ 
ening phase (Section [7. lb . on the post-maximum decline rate and 
the width-luminosity relation l l7.2b . on the late-time evolution d7.3b . 
and on the relative homogeneity of SN la properties in the NIR 
(El. Finally, we comment on the impact of the IGE abundance on 
the near-UV flux level (El- 


7.1 Early brightening phase and high-velocity features 

The early rate of increase in bolometric luminosity is largely insen¬ 
sitive to the ®®Ni mass, as seen through the < 3 per cent variation 
in the power-law exponent a in Tabled The same is true in individ¬ 
ual bands, although it is larger in UBVR and is hardly noticeable 
in I and the NIR bands (Fig. d. This suggests that the early-time 
luminosity, if powered by ®®Ni alone, is more sensitive to mixing 
than slight variations in total ®®Ni mass (see D14a). Indeed, mixing 
can increase the ®®Ni mass fraction by several orders of magnitude 
in the outermost layers (see also Sect. d where an increase in ®®Ni 
mass only results in a modest variation of the outer extent of the 
®®Ni distribution (cf. n(®®Ni) column in Tabled- 

Spectroscopically, the effect is most clearly seen in the UV 
range (blueward of ~3500 A; see Fig. Qo), through variations in 
the absorption by once-ionized IGEs (Sc II, Ti II, Cr II, Ee II, Co II, 
and Nill). The relative strength of the high-velocity components 
to the Call H&K and 8500 A triplet absorptions also varies sig¬ 
nificantly amongst our three models. They are weaker in our model 
with more ®®Ni (DDCIO) due to the higher Ca ionization at 25000- 
35000 km s“ ^. Conversely, the higher excitation in the outer ejecta 
leads to enhanced absorption in the blue wings of several lines, 
including the Eell absorption complex around 4800 A. The in¬ 
crease in ®®Ni mass is not sufficient to match the large Si II 6355 A 
blueshift at early times, whic h might indicate a greater outward 
mixing of ®®Ni (D14a; see also lBaron et alj201^ . 

7.2 Post-maximum decline and the width-luminosity relation 

The bolometric post-maximum decline rate (parametrized via the 
AMis parameter in Tabled changes by less than 0.1 mag from 
model DDCIO to DDC17. Part of this degeneracy with ®®Ni mass 
stems from the increase in the 7 -ray escape fraction for models with 
more ®®Ni. By 20 days past explosion, ~10 per cent of the decay 
energy directly leaks out of the ejecta in model DDCIO, cf. only 
~5 per cent in DDC17. 

However, an increase in the ®®Ni mass results in bluer colours 
at any given time and a more gradual shift to redder colours (Eig. H, 
which affects the post-maximum decline rate in the individual 
bands. In the B-band, the AM 15 parameter varies by ~ 0.4 mag 
(0.92, 1.03, and 1.27 mag in order of decreasing ®®Ni mass), or four 
times the bolometric variation. Despite the relatively narrow range 
in peak Ms, we can compute a width-luminosity relation (WLR) 
for these three models: Mb = 1 . 12 [Ami 5 (B) —1.1] —19.11. This 
WLR is consi s tent w ith the observed one for broad-lined SNe la in 
iBlondin et alj ( l2012h . 

Spectroscopically, the most notable impact of ®®Ni mass 
(other than the absolute flux difference) is in the Cr Il-dominated 
region around 5300 A. As noted in Section [O] this region is par¬ 
ticularly sensitive to the Cr ionization level. 

7.3 Late-time evolution 

The late-time bolometric magnitude decline is the same for our 
three models (^0.028 magday“^), i.e., the rate of change of the 
7 -ray escape fraction is nearly identical. This is a natural conse¬ 
quence of the similar density structure and ®®Ni distribution of the 
inner ejecta in our three Chandrasekhar-mass models. 

The spectra are remarkably similar at these times, with minor 
differences in SED confined to <6000 A, owing to the more uni¬ 
form ionization state of IGEs. Variations in Cr II line opacity con¬ 
tinue to affect the flux around 5300 A. Although not apparent in 
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Wavelength [A] 



Wavelength [A] 


Figure 10. Illustration of the impact of a ±0.1 Mq variation in ®®Ni mass (coloured curves) with respect to our reference model DDC15 (black) on the optical 
(left) and NIR (right) spectra at four selected times between —10 and +60 days from bolometric maximum. The optical spectra have been normalized to the 
same mean flux in the range 5000-6500 A to better visualize relative differences in the SED. The NIR spectra on the other hand ai'e not normalized (but the 
flux has been scaled by for better visibility), such that differences between the models coiTespond to real variations in absolute flux. The tickmarks on the 
ordinate give the zero-flux level. Notice the lower flux level in the near-UV relative to the optical in our model with less ®®Ni. 


Fig-Hnidue to the flux normalization, the [Co III] 5888 A emission 
strength correlates with ®®Ni mass. This correlation is also present 
in observations of S Ne la and supports th e association of this fea¬ 
ture with cobalt (see lPessart et al.|[2014bl . and references therein). 

7.4 NIR homogeneity and secondary maxima 

An increase in ®®Ni mass causes an earlier hardening of the SED 
and a delayed cooling of the ejecta. The result is an earlier decline 
(and hence an earlier first maximum) in the NIR bands (Tabled- 

The associated delayed III—l-II recombination of IGEs has a 
stronger impact in the bluest bands, resulting in a larger magni¬ 
tude variation amongst our three models compared to redder op¬ 
tical and NIR bands (Fig. [TT] see also B13, Appendix C). The 
variation in the I, J, and Ks magnitudes is less than 0.1 mag 
(and occasionally less than 0.05 mag) between the epoch of first 
maximum in those bands up until ~10 days past bolometric max¬ 
imum. Such a small intrinsic dispersion (also reflected in the uni¬ 
formity of their spectra; see Fig. confirms the potential of NIR 
photometry of SNe la for accurate d istance determinations (e.g., 
ikrisciunas. Phillips & Suntzeflll2004h . 

The same effect delays the t ime of NIR secondary max¬ 
ima for larger ®®Ni mass (see also lKas^l2006h . but leaves their 
magnitude unchanged (Tabled, in agreement with observations 



Figure 11. Evolution of the magnitude variation in various bands amongst 
our three models, spanning ~0.2 Mq in ®®Ni mass. 


teiscardi et al.ll2012l . their Fig. 8). The result is a larger magnitude 
variation in the NIR range at 15-20 d from bolometric maximum, 
even exceeding that in the V and R bands. 
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7.5 IGE abundance and near-UV flux level 

The ±0.1 Mq variation in ®®Ni mass represents a non-negligible 
variation in the relative abundances of IGEs and IMEi0 (see Ta¬ 
ble [TJ, which often result in a modest impact on the predicted 
observables, illustrating the degeneracy of SN la properties and 
the difficulties associated with abundance determinations in these 
events (see next section). 

Part of this degeneracy results from variations in ionization 
whose impact often supersedes those associated with variations in 
abundance in determining the ejecta opacity. In the near-UV range, 
where most of the opacity is due to IGEs, it is our model with less 
®®Ni (and hence less IGEs) that has the lowest near-UV flux level 
relative to the optical up until maximum light (Fig. □oj, due to the 
lower ionization level of Ni/Co/Fe. Our models thus show that an 
increase in the IGE abundance i n SNe la does not n ecessarily result 
in more opaque ejecta (see also lSauer et ahll^OOSl) . which compli¬ 
cates the use of near-UV spectra to con strain the progenitor metal- 
licity (see, e.g.. lFolev & Kirshn^l2013h . 

The larger sensitivity of the (7-band flux to the mass of ®®Ni 
results in a < 1 mag variation in this band amongst our three models 
around maximum light (Fig. EJ, and a larger variation mU — V 
compared to o ther colour indic es (Fig. EJ, in line with observational 
findings (e.g. Ijha et al.ll200^ . This large variation could signifi¬ 
cantly affect reddening determinations of SNe la. 


8 COMPARISON TO PREVIOUS WORK 

In-dep th st udies of SN 20Q2bo have been presented bv iBenetti et alj 
( l2004h and IStehle et alj j2005l) . The latter study attempts to infer 
the chemical stratification of the SN 2002bo ejecta based on steady- 
state spectroscopic modelling under the assumption of a sharp pho¬ 
tosphere emitting a pure thermal continuum. Stehle et al. assume 
an ejecta structure (the density profile of W7 shown in Fig. IS and 
vary the abundances in order to reproduce the observed spectra at 
selected epochs. Their abundance profiles are shown in Fig. [m 
For comparison, we overplot the abundance profiles for our model 
DDC15 and those of W7. The elemental yields are reported in Ta¬ 
ble |4] We discuss in turn the inferred chemical stratification and 
the validity of a diffusive inner boundary, which is a fundamental 
ansatz of the method applied bv istehle et al.l l l2005h . 

8.1 Chemical stratification 

The most strikin g difference of our results with those of 
I Stehle et"^ ( l2005h is that they require a large overabundance of 
IGEs and IMEs in the outer layers (> 15000 kms“^) in order to 
reproduce the broad lines of SN 2002bo, including the blue wings 
of Fe II 4800 A at early times. Although our models possess ~ 1000 
times less Ca (by mass) beyond ~20000 km s“^, the resulting syn¬ 
thetic spectra exhibit high-velocity Call features. Such a huge dif¬ 
ference illustrates the difficulties associated with some abundance 
determinations that are highly sensitive to the ionization balance 
and outer density structure. 

In the intermediate layers (10000-15000 kms“^) the distri- 
bution of IMEs is in good agreement with our model, although 
I Stehle et aP ( l2005h infer a larger Mg abundance (0.08 Mq cf. 

® Since all but < 0.01 Mq of the progenitor WD is burnt in our delayed- 
detonation models, the ±0.1 Mq vaiiation in ®®Ni mass results in a com¬ 
parable vaiiation (of opposite sign) in IME mass. 


Table 4. Comparison of nucleosynthetic yields at fexp ~ 0. 


Species 

DDC15 

[Mq] 

Stehle05 

[Mq] 

W7 

[Mq] 

56Ni 

0.511 

0.520 

0.587 

Fe“ 

0.100 

0.360 

0.145 

Ti-I-Cr 

0.018 

0.003 

0.005 

Ca 

0.045 

0.020 

0.012 

S 

0.197 

0.067 

0.085 

Si 

0.306 

0.220 

0.154 

Mg 

0.011 

0.080 

0.013 

Na 

1.68(-5) 

0.001 

6.32(-5) 

0 

0.105 

0.110 

0.138 

C 

0.003 

<0.002 

0.047 


Notes: Numbers in parenthesis correspond to powers of ten. 

“Stable isotopes only (®‘^Fe, ®®Fe, ®^Fe, and ®®Fe). Since these yields are 
given at fexp ~ 0, the stable iron yield does not include ®®Fe from ®®Ni 
decay. 

~ 0.01 Mq in our model) in order to reproduce the Mg II ab¬ 
sorption feature around 4200 A in the early-time spectra. Past 
maximum light, however, the deep absorption at ~4400 A is 
caused by several Fell lines in our model, not by Mg II, hence 
we can reproduce this feature despite the sharp dr op in Mg abun¬ 
dance below 15000 kms~^. Istehle et^ l l2005l) also argue for 
significant downward mixing of oxygen (Xq > 0.1 above ~ 
8500 kms“^), when we satisfactorily fit the Ol 7773 A feature 
up until shortly after maximum light with no oxygen present below 
~ 10000 kms"A _ _ 

Contrary to our model, IStehle et"^ ( l2005l) infer the presence 
of Na and the absence of Ti/Cr below ~ 10000 kms“^. For Na, 
their abundance constraint is most likely a misidentification of the 
emission feature at ~ 5900 A past maximum light, which is due to 
[Com] 5888 A, not Nal (see lPessart et al.ll2014hl) . As for Ti/Cr, 
several CrII lines affect the flux level around ~5300 A past max¬ 
imum light (Sect. 16.31 , so that chromium must be present in the 
inner ejectaQ 

The ^^Ni distribution i n our model is comparable to that in¬ 
ferred bv lStehle~rt ^ (l2005h in the range 5000-10000 km s but 
drops more rapidly beyond 10000 kms“^. With our distribution we 
reproduce the early rise in bolometric luminosity (Fig. 0. Nonethe¬ 
less, we argued for some additional outward mixing of ®®Ni based 
on the lack of absorption in the blue wings of several features in the 
early-time spectra (Sect. | 6 T}. In fact, the ®®Ni distribution of the 
mixed m odel DDCIO M4 of D 14a closely matches the ®®Ni distri¬ 
bution of istehle et alj (l2005l) . Note that we cannot really comment 
on the distribution below ^3000 kms“^ since we did not evolve 
our model beyond ~ 100 days past explosion, when most of the flux 
is expected to emerge from the inner ~0.1 Mq of ejecta. 

8.2 The validity of a diffusive inner boundary 

When modelling the early-time spectra of SN 2002bo IStehle et"^ 
( l2005h assume there exists a sharp photosphere at a pre-defined 
velocity emitting a thermal continuum, yet the importance of line 

At the time their paper was written, IStehle et al.l i2005ll did not have ac¬ 
cess to any spectra of SN 2002bo between +6 d and +263 d past B-band 
maximum, and thus could provide only weak constraints on the abundance 
distribution in these intermediate ejecta layers. 


© 2015 RAS, MNRAS 000.[Tlt35l 



















































Non-LTE spectra of SN 2002bo 13 



5 10 15 20 25 30 35 5 10 15 20 25 30 35 5 10 15 20 25 30 35 

Velocity [10^ km s“^] 


Figur e 12. Comparison of abundance profiles in our reference model DDC15 (red, solid line) with those derived via spectrosc opic modeling bv IStehle et alj 
i2005l) [black; the data were manually scanned from th eir Fig. 5a] and the W7 model of lNomoto. Thielemann & Yokoil il984) [grey]. The ®®Ni distribution 
is shown at fexp ~ 0. Since the ordinate of Fig. 5a of IStehle et all l2005h does not extend below ~3.1X10““^, we give upper limits as dotted lines plus a 
downward-pointing an'ow. We also show the ®®Ni distribution for our models with ±0.1 Mq of ®®Ni (red dotted lines), and for a mixed version of our model 
with -hO.I Mq of ®®Ni (noted DDC10_M4 in D14a; red dashed line). The abundances for Mg through Fe are solar beyond '--'35000 kms“^ in our model. 

scattering in SNe la leads to strong departures from local thermo¬ 
dynamic equilibrium (LTE) in both the gas (level populations) and 
the radiation field (mean intensity Jv)- In our simulations, we solve 
the non-LTE time-dependent radiative transfer for the whole ejecta, 
with a zero-flux inner boundary condition at the base of the inner¬ 
most ejecta shell. We can thus check how .h, and Bi, deviate from 
each other as function of time since explosion, and as a function 
of depth in the ejecta. In Eig. [T^ we illustrate the departures of 
the mean intensity from a Planckian distribution based on our non- 
LTE time-dependent radiative transfer calculations with CMFGEN 
for model DDC15 around bolometric maximum. 

In the innermost regions, at a Rosseland-mean optical depth 
tross > 100 (v < 1500 kms“^), we have « Bv{T) at 
all UV, optical, and NIR wavelengths. Hence, the non-LTE solver 
recovers LTE conditions, i.e., the radiation and the gas are in 
equilibrium at the same (local) temperature T. At tross ~ 30 
(v « 3500 kms“^), Jv < Bu{T) in the red part of the optical 
and in the NIR, but the radiation is still Planckian below '--'SOOO A, 
where the large number of lines ensures an efficient thermaliza- 
tion. As we progress further out to tross of 10 and 1 (~6000 and 
^11000 kms“^, respectively), Jv drops below B^{T) even in the 
UV, although the SED remains smooth below ~2500 A. This per¬ 


sists out to TRoss ~ 0.1 (v ~ 16000 kms“^), although Bi,{T) is 
then no longer a good representation of J^. 

Thus, already at an optical depth of several tens, injecting pho¬ 
tons with a Planckian distribution would lead to an overestimate of 
the flux in the red part of the spectrum. The situation is even worse 
when only the outermost (optically thin) ejecta are treated and one 
as sumes a sharp black body emitting photosphere at r « 1, as done 
bv IStehle et alj ( l2005h . As a result, they invoke a lower reddening 
for SN 2002bo than that inferred from observational constraints. 
In contrast, our full non-LTE solution treats the entire ejecta and 
yields a satisfactory match to the overall SED at all times with a 
total reddening E{B — V) = 0.41 mag (Figs.[8}[9]l. 


9 DISCUSSION AND CONCLUSION 

We have presented a detailed analysis of the photometric and 
spectroscopic evolution of a ID Chandrasekhar-mass delayed- 
detonation model for the Type la SN 2002bo. Beyond providing 
a suitable model for this event, it provides a reference for under¬ 
standing all SNe la similar to that prototype. 

The good match we obtain with the light curves and spectra 
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5000 10000 15000 20000 25000 


Wavelength [A] 


Figure 13. Illustration of the wavelength dependence of the mean intensity 
Ju (solid lines) and the Planck function (dotted lines) at different depths 
in the ejecta for model DDC15 at —2.1 d from bolometric maximum. From 
top to bottom, the panels correspond to locations from the inner to the outer 
ejecta. For each location, we give the velocity, the Rosseland-mean optical 
depth, and the gas temperature. As the optical depth decreases, Ju deviates 
more and more from (T). 


over the first 100 days since explosion gives some support to 
the standard single-degenerate model for these obje cts. It remains 
to be studied whether WD-WD merger models (e.g.. |Pakmor et all 
l2012h can reproduce events like SN 2002bo as well as delayed- 
detonations of single WDs at the Cha ndrasekhar mass. The ejected 
mass distribution of SNe la found bv IScalzo, Ruiter & SimI j20l4) 
is strongly peaked at ~1.4 Mq, yet reveals a significant fraction 
of sub-Chandrasekhar-mass progenitors, although these are asso¬ 
ciated with the least luminous events. While it is unclear to what 
extent a direct confrontation to observations reveals the successes 


and shortcomings either of the radiative-transfer treatment or of 
the underlying explosion model (or both), a proper assessment re¬ 
quires treating the entire SN ejecta in non-LTE at all times (see 
section lS^ . 

Our assumption of spherical symmetry and our somewhat arti¬ 
ficial treatment of the delayed-detonation mechanism are not detri¬ 
mental to reproducing the radiative properties of a proto-typical 
SN la. The mild disagreement between our model and SN 2002bo 
at the earliest times may be resolved by invoking a greater out¬ 
ward mixing of ®®Ni, but could also indicate some fundamental 
characteristic of the explosion not captured by our artificial, one¬ 
dimensional setup. 

In B13 we noted a general overestimate of the width of the 
characteristic Si II 6355 A line in our models to account for the bulk 
of normal SNe la. While this was originally flagged as a potential 
overestimate of the kinetic energy in our input ID hydrodynamical 
models, it has since then become clear that the observed diversity 
in spectral line widths could reflect profound physical differences 
in the explosion mechanism (in particular, pulsating delayed det¬ 
onations result in a narrower Sill 6355 A line; see D14a). The 
delayed-detonation models discussed in B13 are thus well-suited 
to broad-lined SNe la. 

The observed paucity of C II lines in the early-time spectra of 
broad-lined SNe la is then physical and not related, for example, to 
line overlap. In the pulsating delayed-detonation models of D 14a, 
carbon is present down to ~ 15000 kms“^, close to the spectrum- 
formation region at early times. C II lines are thus predicted in 
these models, along with a narrower Si II 6355 A line profile. The 
preferential association of C II detections with narrow-lined SNe la 
(termed “Core Norma l” by Brmichet^ 2006h is backed up by ex¬ 
tensive data sets (e.g.. |PaiTent et alj201 1 ). 

We show that ±0.1 Mq variations in ®®Ni mass have a modest 
impact on the bolometric and colour evolution of our model. Spec¬ 
troscopically, most of the effect is seen blueward of ~3500 A up 
until maximum light, our model with less ®®Ni displaying a lower 
near-UV flux level despite the lower abundance of IGEs. In the 
NIR, an increase in ®®Ni mass delays the secondary maxima but 
does not affect their magnitude, as found in observational studies. 
More importantly, our models confirm the homogeneous nature of 
SNe la in the NIR, with a < 0.1 mag variation in the I, J, and Ks 
bands around maximum light. 

Allowing for homologous expansion of the ejecta inherited 
from the hydrodynamical modelling of the explosion, combined 
with a self-consistent time-dependent radiative-transfer treatment, 
our models yield a better match to SN 2002bo than previous 
parametrized approaches. The complex nature of the radiative 
transfer problem (line overlap, line saturation, uncertain ionization, 
scattering vs. absorption, non-thermal processes, etc.) does not per¬ 
mit a direct and accurate determination of elemental abundances. 
Instead, a global modeling of the ejecta to predict simultaneously 
the bolometric and multi-band light curves, as well as the optical 
and NIR spectra, allows a better assessment of the adequacy of a 
given model. 
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APPENDIX A: CONVERGENCE 

CMFGEN employs an iterative meth od for solving the statistical and 
radiative equilibrium equations tsee lHillieJ 199(11 . A model is con¬ 
verged when the maximum correction to any of the variables in¬ 
volved (temperature, electron density, and level populations) at all 
depths falls below some pre-defined value set by the EPS_TERM 
parameter (typically < 0.1%). The vast majority of corrections are 
well below this limit. 

We illustrate the convergence properties of our DDC15 model 
in Fig. lAll at four selected time steps and for four different ejecta 
locations. We only show the temperature, since its rate of conver¬ 
gence largely determines the global convergence properties of the 
model. Since small corrections to the temperature can result in po¬ 
tentially large corrections of the level populations, the temperature 
corrections are orders of magnitude less than EPS.TERM by the 
time convergence is achieved (typically in 40-70 iterations). 

After the first iteration, the temperature is set to that obtained 
via a grey calculation, and the fractional change is with respect to 
the converged solution of the previous time step, hence the poten¬ 
tially large correction (~ 10%) associated with the first datum. The 
temperature is then held fixed (i.e. the temperature correction is 
zero) during a series of A and full iteratio ns in order to stabilize th e 
corrections to the level populations (see iHillier & DessarlllTOlih . 
Sharp va riations in the temperature correction result from Ng accel¬ 
erations ( lNg|ll974h . used periodically to increase the convergence 
rate. Various diagnostics, such as electron energy balance, global 
energy conservation, and ionization equilibrium, are available as 
independent checks of the convergence. 
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Figure Al. Illustration of the temperature convergence at four selected time steps in our model sequence (increasing from top to bottom). At each time, 
we show the fractional change in temperature |AT|/T = 100 x |Ti_i — Ti\/Ti—i versus iteration number i, at four different ejecta locations: In the 
inner, ^®Ni-deficient region {v 2000 kms~^); in the ^®Ni-rich layers {v ^ 7000 kms~^); in the IME-rich layers (v ^ 12000 kms“^); in the outer 
C/O-dominated layers (u 30000 km s“^). The symbols indicate the sign of the correction (positive for filled upwai'd-pointing triangles; negative for open 
downward-pointing triangles). The dotted line in each plot corresponds to the value of the convergence criterion, noted EPS-TERM in CMFGEN. It is set to 
0.04% at all time steps shown here, except at fexp = 40.11 d where it is set to 0.07%. 
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APPENDIX B: BOLOMETRIC LUMINOSITY AND 
ABSOLUTE MAGNITUDES 

Table HD lists the bolometric and pseudo-bolometric {U Ks) 
luminosities, as well as the optical-NIR absolute magnitudes of 
model DDC15 as a function of time since explosion. 
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Table Bl. Bolometric luminosity, pseudo-bolometric {U Kg) luminosity, and optical-NIR absolute magnitudes for model DDC15. 


^exp 

[d] 

^bol 
[erg s“^] 

[erg s“i] 

U 

[mag] 

B 

[mag] 

V 

[mag] 

R 

[mag] 

I 

[mag] 

J 

[mag] 

H 

[mag] 

Ks 

[mag] 

0.98 

6.02(39) 

4.98(39) 

-5.64 

-7.79 

-9.85 

-11.26 

-11.22 

-12.43 

-12.47 

-12.65 

1.07 

8.71(39) 

7.37(39) 

-5.68 

-8.40 

-10.45 

-11.77 

-11.67 

-12.74 

-12.80 

-12.93 

1.18 

1.27(40) 

1.10(40) 

-6.11 

-9.00 

-11.05 

-12.27 

-12.14 

-13.05 

-13.12 

-13.22 

1.30 

1.84(40) 

1.63(40) 

-6.65 

-9.61 

-11.62 

-12.74 

-12.61 

-13.36 

-13.44 

-13.52 

1.43 

2.62(40) 

2.37(40) 

-7.38 

-10.26 

-12.16 

-13.14 

-13.01 

-13.65 

-13.75 

-13.82 

1.57 

3.69(40) 

3.37(40) 

-8.28 

-10.91 

-12.68 

-13.52 

-13.37 

-13.95 

-14.05 

-14.11 

1.73 

5.15(40) 

4.72(40) 

-9.24 

-11.54 

-13.16 

-13.87 

-13.71 

-14.24 

-14.33 

-14.39 

1.90 

7.06(40) 

6.48(40) 

-10.13 

-12.11 

-13.60 

-14.18 

-14.01 

-14.51 

-14.60 

-14.65 

2.09 

9.65(40) 

8.85(40) 

-10.93 

-12.63 

-14.03 

-14.48 

-14.31 

-14.78 

-14.87 

-14.93 

2.30 

1.31(41) 

1.20(41) 

-11.61 

-13.11 

-14.42 

-14.76 

-14.61 

-15.06 

-15.14 

-15.20 

2.53 

1.76(41) 

1.60(41) 

-12.23 

-13.55 

-14.79 

-15.04 

-14.89 

-15.33 

-15.40 

-15.46 

2.78 

2.33(41) 

2.13(41) 

-12.84 

-13.98 

-15.11 

-15.32 

-15.16 

-15.59 

-15.64 

-15.72 

3.06 

3.11(41) 

2.84(41) 

-13.43 

-14.40 

-15.43 

-15.60 

-15.43 

-15.86 

-15.89 

-15.98 

3.37 

4.17(41) 

3.83(41) 

-14.01 

-14.84 

-15.74 

-15.89 

-15.70 

-16.13 

-16.15 

-16.24 

3.70 

5.67(41) 

5.20(41) 

-14.57 

-15.29 

-16.06 

-16.20 

-15.99 

-16.42 

-16.42 

-16.5! 

4.07 

7.63(41) 

7.05(41) 

-15.07 

-15.73 

-16.35 

-16.49 

-16.26 

-16.68 

-16.68 

-16.76 

4.48 

1.02(42) 

9.50(41) 

-15.54 

-16.17 

-16.64 

-16.78 

-16.53 

-16.95 

-16.93 

-16.99 

4.93 

1.36(42) 

1.27(42) 

-15.99 

-16.58 

-16.91 

-17.06 

-16.79 

-17.19 

-17.16 

-17.21 

5.42 

1.78(42) 

1.67(42) 

-16.43 

-16.96 

-17.16 

-17.33 

-17.04 

-17.43 

-17.37 

-17.39 

5.96 

2.32(42) 

2.19(42) 

-16.87 

-17.31 

-17.41 

-17.58 

-17.30 

-17.64 

-17.56 

-17.56 

6.56 

2.98(42) 

2.83(42) 

-17.31 

-17.63 

-17.64 

-17.82 

-17.58 

-17.83 

-17.73 

-17.70 

7.22 

3.78(42) 

3.60(42) 

-17.72 

-17.92 

-17.85 

-18.04 

-17.84 

-17.99 

-17.87 

-17.83 

7.94 

4.70(42) 

4.48(42) 

-18.09 

-18.17 

-18.06 

-18.25 

-18.07 

-18.14 

-17.98 

-17.93 

8.73 

5.70(42) 

5.42(42) 

-18.39 

-18.39 

-18.25 

-18.43 

-18.25 

-18.26 

-18.06 

-18.00 

9.60 

6.74(42) 

6.40(42) 

-18.65 

-18.59 

-18.43 

-18.59 

-18.40 

-18.36 

-18.10 

-18.05 

10.56 

7.80(42) 

7.40(42) 

-18.88 

-18.76 

-18.60 

-18.73 

-18.52 

-18.43 

-18.12 

-18.06 

11.62 

8.78(42) 

8.34(42) 

-19.04 

-18.89 

-18.75 

-18.86 

-18.62 

-18.51 

-18.13 

-18.07 

12.78 

9.63(42) 

9.22(42) 

-19.15 

-19.01 

-18.88 

-18.97 

-18.67 

-18.53 

-18.10 

-18.04 

14.06 

1.05(43) 

9.90(42) 

-19.31 

-19.10 

-18.97 

-19.05 

-18.72 

-18.53 

-18.04 

-17.98 

15.47 

1.09(43) 

1.05(43) 

-19.36 

-19.14 

-19.05 

-19.12 

-18.74 

-18.50 

-17.99 

-17.94 

17.02 

1.16(43) 

1.08(43) 

-19.42 

-19.14 

-19.17 

-19.17 

-18.83 

-18.45 

-17.98 

-18.07 

18.72 

1.13(43) 

1.07(43) 

-19.31 

-19.12 

-19.21 

-19.19 

-18.83 

-18.35 

-17.94 

-18.05 

20.59 

1.08(43) 

1.03(43) 

-19.18 

-19.07 

-19.23 

-19.19 

-18.81 

-18.24 

-17.95 

-18.06 

22.65 

1.00(43) 

9.70(42) 

-19.02 

-18.99 

-19.21 

-19.16 

-18.77 

-18.09 

-18.00 

-18.08 

24.91 

9.01(42) 

8.72(42) 

-18.82 

-18.85 

-19.16 

-19.08 

-18.70 

-17.81 

-18.06 

-18.1! 

27.40 

7.76(42) 

7.51(42) 

-18.57 

-18.65 

-19.05 

-18.95 

-18.59 

-17.33 

-18.09 

-18.13 

30.14 

6.56(42) 

6.36(42) 

-18.29 

-18.39 

-18.89 

-18.79 

-18.52 

-17.05 

-18.18 

-18.17 

33.15 

5.74(42) 

5.54(42) 

-17.83 

-17.97 

-18.70 

-18.70 

-18.72 

-17.34 

-18.62 

-18.54 

36.46 

5.14(42) 

4.93(42) 

-17.43 

-17.58 

-18.47 

-18.61 

-18.85 

-17.61 

-18.80 

-18.70 

40.11 

4.47(42) 

4.15(42) 

-17.05 

-17.18 

-18.17 

-18.45 

-18.83 

-17.92 

-18.80 

-18.70 

44.12 

3.51(42) 

3.28(42) 

-16.75 

-16.86 

-17.85 

-18.17 

-18.58 

-17.76 

-18.53 

-18.41 

48.53 

2.74(42) 

2.56(42) 

-16.50 

-16.61 

-17.59 

-17.91 

-18.30 

-17.35 

-18.22 

-18.05 

53.38 

2.20(42) 

2.05(42) 

-16.31 

-16.42 

-17.38 

-17.69 

-18.03 

-16.92 

-17.92 

-17.72 

58.72 

1.80(42) 

1.67(42) 

-16.15 

-16.27 

-17.20 

-17.49 

-17.76 

-16.50 

-17.64 

-17.43 

64.59 

1.48(42) 

1.37(42) 

-16.00 

-16.14 

-17.03 

-17.30 

-17.49 

-16.08 

-17.37 

-17.15 

71.05 

1.22(42) 

1.13(42) 

-15.84 

-16.02 

-16.85 

-17.10 

-17.21 

-15.65 

-17.09 

-16.87 

78.16 

1.00(42) 

9.29(41) 

-15.68 

-15.90 

— 16.66 

-16.89 

-16.94 

-15.20 

-16.80 

-16.59 

86.00 

8.15(41) 

7.56(41) 

-15.49 

-15.77 

-16.45 

-16.65 

— 16.66 

-14.74 

-16.48 

-16.29 

94.60 

6.57(41) 

6.10(41) 

-15.29 

-15.63 

-16.23 

-16.39 

-16.38 

-14.26 

-16.12 

-15.98 

104.10 

5.26(41) 

4.86(41) 

-15.05 

-15.48 

-15.98 

-16.10 

-16.10 

-13.8! 

-15.73 

-15.66 


Note: Numbers in p arentheses correspond to powers of ten. UBVRI magnitudes are based on the passbands of iBesse iiC^; JHKs magnitudes are in 
the 2MASS system jCohen. Wheaton & Megeathll2003h . 
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APPENDIX C: STRONG LINES 


Table Cl. Strong lines in the optical and NIR at 4.1 d past explosion. 


^exp 

= 4.1d(Aiboi = 

-13.5d) 

Ion 

3000-10000 A 
-^air [A] 

1 .0-2.5 fim 
Aair [/rm] 

Oi 

9266.01 

1.316 

Mg II 

7896.37 

1.009 

AllI 

7042.08 

1.091 

1.863 

2.137 

2.404 

1.008 

Si II 

5978.93 

1.175 

Sim 

6347.11 

9412.63 

4828.95 

1.370 

1.445 

1.691 

1.947 

1.250 

SlI 

5739.73 

4162.67 

1.523 

Sill 

6755.65 


Am 

7171.58 

4348.06 


Call 

3933.66t 

1.184 

Ti II 

8542.09 

3088.04 


Cm 

3234.52 

3349.41 

3759.30 

3124.97 


Fell 

3002.65 


Felll 

4419.60 


Coll 

5156.11 

3501.72 


Nill 

3513.99 



den and semi-forbidden transitions are noted using the appropriate 
brackets around the line wavelength (in air). We highlight transi¬ 
tions connected to the ground state with a “f ” symbol. These tables 
are meant to serve as a guide for line identification in spectra of 
SNe la similar to SN 2002bo. However, we caution the reader that 
the Sobolev EWs are only approximate and neglect important ef¬ 
fects such as line overlap. Moreover, lines with a P-Cygni profile 
morphology may be important contributors to the resulting spec¬ 
trum yet yield a small absolute EW. This explains for instance why 
the strong Call K line (3934 A) does not show up in Table C2 
(fexp = 5.4 d) where Fig. D1 shows the strong impact of the 
Call H&K feature on the optical spectra at these times, and the 
Call K line indeed appears in Tables Cl (fexp = 4.1 d) and C3 
(fexp = 9.6 d). 


In tables IcTl - IClSI we report the strongest lines in the optical (3000- 
10000 A) and NIR (1-2.5 fim) ranges at a given time between 4.1 d 
and 104.1 d past explosion (—13.5 d and -1-86.5 d from bolometric 
maximum). A line is considered strong when its absolute Sobolev 
equivalent width (EW) exceeds 10 per cent of the largest absolute 
EW within ±10^ kms“^ of the line’s central wavelength. Forbid- 


(c) 2015 RAS, MNRAS 00Q.[Tll35l 








Non-LTE spectra of SN 2002bo 21 


Table C2. Strong lines in the optical and NIR at 5.4 d past explosion. 


Table C3. Strong lines in the optical and NIR at 9.6 d past explosion. 


^exp 

= 5.4d(Atboi = 

-12.Id) 

Ion 

3000-10000 A 

-^air [^] 

1 .0-2.5 fim 
Aair [Aim] 

Oi 

7771.94 

1.316 

Mg II 

7896.37 

1.091 

AlII 

7042.08 

1.863 

2.137 

2.404 

1.008 

Sill 

5978.93 

1.175 

Si III 

6347.11 

9412.63 

4552.62 

1.370 

1.445 

1.691 

1.250 

SlI 

4828.95 

5739.73 

7466.32 

8262.57 

9799.91 

5453.85 

1.502 

1.318 

Sill 

3497.28 

1.032 

Am 

3838.27 

4253.50 

5219.31 

6755.65 

7171.58 

4806.02 


Call 

8542.09 

1.856 

Sell 

3572.53 


Till 

3088.04 


Cm 

3234.52 

3124.97 


Fell 

3002.65 


Felll 

5156.11 


Nill 

3087.07 



^exp 

= 9.6d (Atbol = 

-8.0 d) 

Ion 

3000-10000 A 
Aair [A] 

1.0-2. 5 fim 
Aair [Aim] 

Oi 

7771.94 

1.316 

Mg II 

9218.25 

1.091 

Si II 

6347.11 

1.370 

Sim 

4552.62 

1.691 

2.199 

1.134 

Si IV 

4828.95 

5739.73 

7466.32 

8262.57 

9323.89 

9799.91 

6701.21 

1.250 

1.502 

2.093 

1.190 

SlI 

7720.37 

8957.25 

5453.85 

1.855 

1.985 

2.426 

1.318 

Sill 

3234.01 

1.032 

SIV 

3497.28 

4253.50 

5219.31 

6140.66 

1.421 

Am 

7003.27 

4806.02 


Ami 

7988.48 

1.312 

Call 

3933.66t 

1.856 

Till 

8542.09 

3349.41 


Cm 

3124.97 


Fell 

5169.03 


Felll 

3013.13 


Coll 

3501.72 


Com 

3010.88 


Nill 

3513.99 
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Table C4. Strong lines in the optical and NIR at 12.8 d past explosion. 


Table C5. Strong lines in the optical and NIR at 15.5 d past explosion. 


^exp 

= 12.8d(Atboi 

= -4.8 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [^m] 

Oi 

7771.94 


Mg II 

9218.25 

1.091 

Sill 

5978.93 

1.370 


6347.11 

1.691 

2.199 

Sim 

4552.62 

1.134 


5739.73 

1.250 


7466.32 

1.502 


8262.57 

9323.89 

2.093 

Si IV 

6701.21 

1.190 


7720.37 

1.855 


8957.25 

1.985 

2.426 

SlI 

4815.55 

5453.85 

1.362 

Sill 

3497.28 

1.032 


4253.50 

5219.31 

1.421 

SiV 

7003.27 


Am 

4806.02 


Arm 

[7135.79]t 

1.312 

Call 

3933.66t 

8542.09 

1.856 

Till 

3349.41 


Cm 

3124.97 


Fell 

5169.03 


Felll 

3013.13 


Coll 

3501.72 

2.221 

Com 

3010.88 

3305.27 

[2.281] 

Co IV 

8033.21 


Nill 

3513.99 



^exp 

= 15.5d(Atboi 

= -2.1 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [Aim] 

Oi 

7771.94 


Mg II 

9218.25 

1.091 

Sill 

6347.11 

1.691 

2.199 

Sim 

4552.62 

1.134 


5739.73 

1.250 


7466.32 

9323.89 

1.502 

Si IV 

6701.21 

1.190 


7720.37 

1.855 

2.426 

SlI 

5453.85 

1.362 

Sill 

3497.28 

1.032 


4253.50 

5219.31 

1.421 

SiV 

7003.27 


Am 

4806.02 

1.047 

Ami 

[7135.79]t 

1.440 

Call 

3933.66t 

8542.09 

1.856 

Till 

3349.41 


Cm 

3124.97 


Fell 

5169.03 


Felll 

3013.13 


FelV 


1.868 

Coll 

3501.72 

2.221 

Com 

3010.88 

[2.002] 


3305.27 

[2.097] 


4740.66 

[2.281] 

Co IV 

8033.21 


Nill 

3513.99 


Nilll 

4874.95 
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Table C6. Strong lines in the optical and NIR at 17.0 d past explosion. 


Table C7. Strong lines in the optical and NIR at 22.6 d past explosion. 


^exp 

= 17.0d(Atboi 

= -0.5 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [^m] 

Oi 

7771.94 


Mg II 

7896.37 

1.091 

Sill 

6347.11 

1.175 



1.370 



1.691 



2.199 

Sim 

4552.62 

1.134 


5739.73 

7466.32 

9323.89 

1.250 

SlI 

5453.85 

1.362 

Sill 

3497.28 

1.032 


4253.50 

1.421 


5219.31 

6755.65 

2.374 

Am 

4806.02 

1.047 

Arm 

[7135.79]t 

1.312 

Call 

3933.66t 

8542.09 

1.885 

Cm 

3124.97 


Fell 

5169.03 


Felll 

5156.11 

[2.218] 

Coll 

3501.71 

1.777 

Com 

3010.88 

[1.548] 


3305.27 

[2.002] 


4740.66 

[2.097] 


5544.79 

[2.281] 


6078.55 

[2.473] 

Nill 

3513.99 


Nilll 

[7889.93] 

1.001 


^exp 

= 22.6d(Atbol 

= +5.1 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [Aim] 

Oi 

7771.94 


Mg II 

9218.25 

1.091 

Sill 

6347.11 

1.175 



1.445 



1.691 

Sim 

5739.73 

1.134 


9323.89 

1.250 

SlI 

5453.85 

[1.032] 

Sill 

4253.50 

1.087 

Am 

6643.69 


Ami 

[7135.79]t 


Call 

3933.66t 

8542.09 

1.184 

Till 

3349.41 


Cm 

3132.05 


Fell 

3002.65 

1.086 

Felll 

4419.60 

5156.11 

1.365 

FelV 

[7925.02] 

[1.076] 

Coll 

3621.20 

1.606 

Com 

3151.40] 

[1.548] 


3259.60 

[2.002] 


3782.18 

[2.097] 


4740.66 

[2.281] 


6994.09 

[2.473] 

Co IV 

[7902.03] 

[1.871] 

[2.382] 

Nill 

3513.99 


Nilll 

[7889.93] 

1.001 
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Table C8. Strong lines in the optical and NIR at 27.4 d past explosion. 


^exp 

= 27.4d(Atbol 

= +9.8 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [^m] 

Oi 

7771.94 


Mg II 

9218.25 

1.091 

Sill 

4130.90 

1.175 


6347.11 

1.370 



1.445 

Sim 

4552.62 

1.134 

SlI 

5453.85 

[1.032] 

Sill 

[9530.62] 


Am 

6643.69 


Arm 

[7135.79]t 


Call 

3933.66t 

1.885 


8542.09 


Till 

3349.41 


Cm 

3132.05 


Fell 

5169.03 

1.113 

Felll 

4431.02 

1.365 


7265.19 


Coll 

3501.71 

1.606 


3621.20 

1.777 



2.220 



2.361 

Com 

5667.05 

[1.272] 


[5888.48]t 

[1.548] 



[2.002] 



[2.097] 



[2.473] 

Co IV 


[1.871] 

Nill 

3513.99 

1.878 

Nilll 

[7889.93] 

1.001 


Table C9. Strong lines in the optical and NIR at 33.2 d past explosion. 


^exp 

= 33.2d(Atbol 

= +15.6d) 

Ion 

3000-10000 A 
Aair [A] 

1 .0-2.5 fim 
Aair [+m] 

Si II 

6347.11 

1.175 

So 

5647.02 

[1.032] 

Sm 

[9530.62] 


Arm 

[7135.79]t 


Call 

3968.47t 

1.184 

Till 

8542.09 

3349.4! 


Cm 

3132.05 


Fell 

5169.03 

1.050 

Felll 

7462.41 

7711.72 

9997.58 

4419.60 

[1.061] 

Coll 

3621.20 

1.136 

Com 

4160.66 

4660.63 

[5888.48]t 

1.183 

1.350 

1.406 

1.606 

1.777 

1.886 

2.134 

2.220 

2.361 

2.460 

[1.272] 

Co IV 


[2.002] 

[1.871] 

Nill 

3513.99 

1.878 

Nilll 

[7889.93] 
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Table CIO. Strong lines in the optical and NIR at 44.1 d past explosion. 


Table Cll. Strong lines in the optical and NIR at 53.4 d past explosion. 


^exp 

= 44.1d(Atboi 

= +26.6 d) 

Ion 

3000-10000 A 
Aair [A] 

1 .0-2.5 fim 
Aair [^m] 

Si n 

6347.11 


Sll 

[4068.60]t 

[1.032] 

Sill 

[9530.62] 


Arm 

[7135.79]t 


Call 

3968.47t 

1.184 

Till 

8542.09 

8662.14 

3372.80 


Cm 

3132.05 


Fell 

4351.77 

1.050 

Fell! 

4549.47 

5169.03 

7462.41 

7711.72 

9997.58 

[4658.05]t 

1.261 

Coll 

3621.20 

1.183 

Com 

9418.11 

[5888.48]t 

1.406 

1.503 

1.606 

1.777 

1.886 

2.134 

2.220 

2.361 

2.460 

[2.002] 

Nill 

8096.75 

1.312 


^exp 

= 53.4d(Atbol 

= +35.8 d) 


3000-10000 A 

1 .0-2.5 fim 

Ion 

Aair [A] 

Aair [+m] 

Si ll 

6347.11 


Sll 

[4068.60]t 

[1.032] 

Sm 

[9530.62] 


Arm 

[7135.79]t 


Call 

3179.33 

3968.47t 

8542.09 

1.184 

Till 

3372.80 

3759.30 

1.463 

Fell 

4233.17 

1.050 


4549.47 

1.261 


5169.03 

1.370 


7462.41 

7711.72 

9997.58 

1.453 

Felll 

[4658.05]t 


Coll 

3621.20 

1.183 



1.606 



1.777 



1.886 



2.134 



2.220 



2.361 



2.460 

Com 

[5888.48]t 

[2.002] 

Nill 


[1.939] 


© 2015 RAS, MNRAS OOO.rTHBSl 










26 S. Blondin et al. 


Table C12. Strong lines in the optical and NIR at 64.6 d past explosion. 


Table C13. Strong lines in the optical and NIR at 71.1 d past explosion. 


^exp 

= 64.6d(Atbol 

= +47.0 d) 

^exp 

= 71.1d(Atbol 

= +53.5 d) 

Ion 

3000-10000 A 
Aair [A] 

1.0-2.5 fim 

Aair [+m] 

Ion 

3000-10000 A 
Aair [A] 

1.0-2.5 fim 
Aair [+m] 

Si n 

6347.11 


Si II 

6347.11 


Sll 

[4068.60]t 

[1.032] 

Sll 

[4068.60]t 

[1.032] 

Sill 

[9068.62] 


Sm 

[9068.62] 



[9530.62] 



[9530.62] 


Arm 

[7135.79]t 


Arm 

[7135.79]t 


Call 

3933.66t 

1.184 

Call 

3933.66t 

1.184 


8542.09 



8542.09 


Tin 

3759.30 

1.463 

Till 

3759.30 

1.463 

Cm 

3132.05 

1.880 

Cm 

3132.05 

1.503 


3368.04 



3368.04 

1.880 


3408.76] 



3408.76] 


Mnll 

3441.99 


Mnll 

3441.99 


Fell 

4233.17 

1.050 

Fell 

4233.17 

1.050 


4549.47 

1.261 


4549.47 

1.261 


5169.03 

1.370 


5169.03 

1.370 


6247.56 

1.453 


6247.56 

1.687 


6456.38 



6456.38 

1.840 


7462.41 



7462.41 



7711.72 



7711.72 



9997.58 



9997.58 


Felll 

[4658.05]t 

[2.145] 

Felll 

[4658.05]t 

[2.145] 

Coll 

3621.20 

1.183 

Coll 

3621.20 

1.183 



1.606 



1.606 



1.777 



1.777 



2.134 



2.134 



2.220 



2.220 



2.361 



2.361 



2.460 



2.460 

Com 

[5888.48]t 

[2.002] 

Com 

[5888.48]t 

[2.002] 

Co IV 

[3936.51]t 

[1.871] 

Co IV 

[3936.51]t 

[1.871] 

Nill 

[7377.83]t 

[1.939] 

Nill 

[7377.83]t 

[1.939] 

Nilll 

[7889.93] 


Nilll 

[7889.93] 
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Table C14. Strong lines in the optical and NIR at 86.0 d past explosion. 


Table C15. Strong lines in the optical and NIR at 104.1 d past explosion. 


^exp 

= 86.0d(Atboi 

= +68.4 d) 


3000-10000 A 

1.0-2.5 fim 

Ion 

Aair [A] 

Aair [+m] 

Si n 

6347.11 


Sll 

[4068.60]t 

[1.032] 

Sill 

[9068.62] 

[9530.62] 


Arm 

[7135.79]t 


Call 

3933.66t 

8542.09 

1.184 

Tin 

3759.30 

1.463 

Cm 

3132.05 

3368.04 

3408.76] 

4558.65 

1.454 

Mnll 

3441.99 


Fel 

3581.19 

3734.86 

1.197 

Fell 

4233.17 

1.050 


5169.03 

[1.257]t 


6247.56 

1.370 


6456.38 

1.687 


7462.41 

7711.72 

9997.58 

1.840 

Felll 

[4658.05]t 

[2.145] 

Coll 

3501.71 

1.183 



1.606 



1.777 



2.220 



2.361 



2.460 

Com 

[5888.48]t 

[1.548] 


[6853.53] 

[2.002] 

[2.097] 

Co IV 

[3936.51]t 

[4494.29] 

[1.871] 

Nill 

[7377.83]t 

[1.939] 

Nilll 

[7889.93] 



^exp 

= 104.1 d(Atboi 

= +86.5 d) 

Ion 

3000-10000 A 
Aair [A] 

1.0-2.5 fim 
Aair [+m] 

So 

[4068.60]t 

[1.032] 

Sm 

[9068.62] 


Arm 

[9530.62] 

[7135.79]t 


Call 

3933.66t 

1.184 

Till 

8542.09 

3759.30 

1.463 

Cm 

4395.03 

3132.05 

1.763] 

Mnll 

3368.04 

3408.76] 

3441.99 


Fel 

3581.19 

1.197 

Fell 

3734.86 

4233.17 

[1.257]t 

Felll 

5169.03 

6247.56 

6456.38 

7462.41 

7711.72 

[4658.05]t 

[1.321] 

1.370 

1.687 

1.840 

[2.145] 

Coll 

3501.71 

[2.218] 

[2.348] 

1.606 

Com 

[5888.48]t 

1.777 

2.460 

[1.548] 

Co IV 

[6853.53] 

[4494.29] 

[2.002] 

[1.871] 

Nill 

[7411.61] 

[1.939] 

Nilll 

[7889.93] 
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APPENDIX D: CONTRIBUTION OF INDIVIDUAL IONS 
TO THE TOTAL OPTICAL AND NEAR-INFRARED FLUX 

Figures lDT4D7l reveal the contribution of individual ions (bottom 
panels) to the full optical (left) and NIR (right) synthetic spectra of 
DDC15 (top panels, red line), compared to SN 2002bo (top panels, 
black line) at all times presented in Figures]^ andOnly ions that 
impact the flux at the > 10 per cent level are shown. 
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Wavelength [A] 



Wavelength [A] 



Wavelength [A] 



Wavelength [A] 


Figure Dl. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at —12.9 d and —11.7 d from from pseudo-bolometric {U —>■ Ks) maximum. 
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Wavelength [A] 



Wavelength [A] 



Wavelength [A] 



Wavelength [A] 


Figure D2. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at —8.5 d (NIR only), —7.7 d (optical only), and —4.1 d from from pseudo-bolometric {U —>■ Ks) maximum. 
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Wavelength [A] 



Wavelength [A] 



Wavelength [A] 



Wavelength [A] 


Figure D3. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at —2.1 d and +0.1 d from from pseudo-bolometric {U Ks) maximum. 
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Wavelength [A] 



Wavelength [A] 



10000 12000 14000 16000 18000 20000 22000 24000 

Wavelength [A] 



Wavelength [A] 


Figure D4. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at +5.9 d, +11.1 d (optical only) and +11.2 d (NIR only) from from pseudo-bolometric {U —)• Ks) maximum. Ion 
spectra marked with a have been scaled down for clarity. 
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Wavelength [A] 



Wavelength [A] 



Wavelength [A] 



Wavelength [A] 


Figure D5. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at +16.2 d, +28.2 d (optical only), and +29.3 d (NIR only) from from pseudo-bolometric {U —>■ Ks) maximum. Ion 
spectra marked with a have been scaled down for clarity. 
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3000 4000 5000 6000 7000 8000 9000 10000 


Wavelength [A] 



3000 4000 5000 6000 7000 8000 9000 10000 


Wavelength [A] 



Wavelength [A] 



Wavelength [A] 


Figure D6. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at +45.1 d, +56.3 d (NIR only), and +56.9 d (optical only) from from pseudo-bolometric {U —>■ Ks) maximum. Ion 
spectra marked with a have been scaled down for clarity. 
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Figure D7. Contribution of individual ions (bottom panels) to the full optical (left) and NIR (right) synthetic spectra of DDC15 (top panels, red line), compared 
to SN 2002bo (top panels, black line) at +72.9 d, +83.6 d (optical only), and +84.6 d (NIR only) from from pseudo-bolometric {U —>■ Ks) maximum. Ion 
spectra marked with a have been scaled down for clarity. 
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